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ABSTRACT 

We present 3.6 and 4.5 /^m Spitzer IRAC imaging over 0.77 square degrees at the Virgo cluster core for 
the purpose of understanding the formation mechanisms of the fow surface brightness intracluster hght 
features. Instrumental and astrophysical backgrounds that are hundreds of times higher than the signal 
were carefully characterized and removed. We examine both intracluster light plumes as well as the outer 
halo of the giant elliptical M87. For two intracluster light plumes, we use optical colors to constrain 
their ages to be greater than 3 & 5 Gyr, respectively. Upper limits on the IRAC fluxes constrain the 
upper limits to the masses, and optical detections constrain the lower limits to the masses. In this first 
measurement of mass of intracluster light plumes we find masses in the range of 5.5 x 10^ — 4.5 x 10^ 
and 2.1 x 10^ — 1.5 x 10^ Mq for the two plumes for which we have coverage. Given their expected short 
lifetimes, and a constant production rate for these types of streams, integrated over Virgo's lifetime, 
they can account for the total ICL content of the cluster implying that we do not need to invoke ICL 
formation mechanisms other than gravitational mechanisms leading to bright plumes. We also examined 
the outer halo of the giant elliptical M87. The color profile from the inner to outer halo of M87 (160 
Kpc) is consistent with either a flat or optically blue gradient, where a blue gradient could be due to 
younger or lower metallicity stars at larger radii. The similarity of the age predicted by both the infrared 
and optical colors (> few Gyr) indicates that the optical measurements are not strongly affected by dust 
extinction. 

Subject headings: galaxies: clusters: individual — galaxies: evolution — galaxies: photometry — 
infrared: galaxies — cosmology: observations 



1. INTRODUCTION 

The basic picture of hierarchical structure formation 
is that large galaxies assemble from smaller pieces. A 
galaxy's environment during assembly is expected to play 
a significant role in its evolution. Indeed, galaxies in 
the dense cores of clusters differ systematically from 
those in the coeval field, in terms of their morphologies, 
stellar populations, and gas co ntent ( e.s. iDressler et all 
1997t iButcher fc Oemled Il978l. 119841: iDesai et all l2007t 
Krick et all 120091) ! 

Many different physical processes have been implicated 
in driving the evolution of cluster galaxies, including both 
gravitational and gas dynamical effects. Gravitational ef- 
fects includesfow galaxy-gala xy interactions and merg- 
ers ()Mihosl[2003 ). harassment (iMoore et al.lll996l Il998li 
the e ffects of the global tidal field (e.g. iBvrd fc ValtonenI 
119901 : iHcnriksen fc By rd 1996) and the effects of clus - 
ter sub-cluster merging IBekki 1 995 IStevens et al.lll999tl . 
Gas dynamical effects include ram pr essure and turbulent 
viscous stripping of the cold ISM ( Gunn fc Gott 1972t 
Abadi et al.l 1199 9: Schu lz fc Struckl boOlt IVollmer et al ' 
20011 : iMurpiiTet al.. .2009^ and starvation (ILarson et al 
1980). 



Within a cluster, gravitational effects can lead to the 
removal of stars from their parent galaxies into the intr- 
acluster medium (ICM). Gas dynamical effects can addi- 
tionally strip gas and dust from galaxies. This stripped 



material may subsequently become sites of star forma- 
tion within the ICM. The properties of these stars, as 
traced by the intracluster light (ICL), can constrain the 
types and f requency of the p h ysical processes at work 
in clusters (|K rick et al. 2006; Krick & Bernstein 2 007 



Gonzalez et al.i 2005,; Fcldnicicr ct_aL 2004; Zibctti et al.l 
20051) by serving as a fossil record of the interaction his- 



tory of the cluster. 

Large stellar streams from galaxy interactions are a 
gravitational effect we know occurs in clusters. Massive 
streams will be visible for of order ~ 1 — 2 times their dy- 
namic al time before dispe rsing into the diffuse ICL popu- 
lation (iRudic k et al.l[2009l ). If the streams can account for 
all of the diffuse ICL built up over time, then there is no 
need to invoke other mechanisms such as gas dynamical 
effects including ram pressure stripping or starvation in 
generating large amounts of ICL. 

Virgo is one of the best places to study these evolution- 
ary processes in action. First, Virgo is the only clu ster 
that is close enough (16.6 Mpc; IShaplev et all l200lD to 
easily achieve many resolution elements per square kilopar- 
sec. Second, many of its g alaxies have clearly be en modi- 
fied by their environment (lAbramson et all 120 111) . Third, 
the cluster exhibits substantial substructure, with infalling 
groups containing tidally interacting galax ies and mergers, 
as se en for the M87 fc M 86 regions (e.g. ISchindler et al.l 
1999; iKennev et alll2008[ ). Finally, Virgo is well-studied 
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at a range of wavelengths (there are at least 677 literature 
references in NED for the Virgo cluster). 

As a central dominant (cD) galaxy of the Virgo cluster, 
M87 is a prototypical giant elliptical. Its optical surface 
brightness profile is poorly fit by a single deV profile; in- 
stead , better fits are achieved through high-n Sersi c pro- 
files (jKormendv et al.l l2009t iJanowiecki et al.l |2010() or a 
combination of single deV profile with an additional outer 
component (Liu et al 2005). Either profile shape implies 
a significant amount of luminosity at large radius, possi- 
bly attributable to an ICL component. M87 is known to 
have a flat to slightly blue c olor gradient from the inside to 
outer regions of the galaxy (|Liu et al.ll2005l : IZeilinger et al.l 
Il993t ). A flat color profile implies a constant age as a func- 
tion of radius. A color gradient, even in the near infrared, 
would imply changing age as a function of radius. Flat 
age profiles cannot rule out monolithic collapse scenarios 
for the formation of this giant elliptical . 

At present, the only cluster- wide information obtained 
for the distribution of ICL in Vir go comes from the opti- 
cal (VkB-hand) ob servations of Mihos et al.l (|2005[) and 
iRudick et al.l (|2010[ ). These data reveal an intricate com- 
plex of intracluster stellar features, including long (> 
lOOkpc) streamers, arcs, tails, bridges, and diffuse blobs 
unassociated with galaxies. These data are clear evi- 
dence for interactions between the infalling populations 
of smaller galaxies with the cluster cD (M87) and the in- 
falling M86 sub-cluster structure. 

In this paper we present a map of the central 0.77 square 
degrees in Virgo using IRAC at 3.6 and 4.5 /im. These 
wavelengths are interesting because they 1) are sensitive 
estimators of stellar mass, including that contributed by 
old stars; 2) will provide a color measurement of the dif- 
fuse features in this cluster; and 3) are largely unaffected 
by the presence of dust. 

This paper is structured in the following manner. In 
fj2| and ij3| we discuss the data including specific reduc- 
tion techniques for low surface brightness (LSB) measure- 
ments with Spitzer. Detection limits are described in 
m Sj5| discusses the results for the intracluster plumes 
and halo of M87. The paper is summarized and conclu- 
sions are drawn in ^ Throughout this paper we use 
Ho = 7Ikm/s/ Mpc and = 0-27 . With this cosmol- 
ogy, we use the Shaplcv et al.l (|2P0I) luminosity distance 
of 16.6 Mpc corresponding to a redshift of 0.0039. 

2. SPITZER WARM IRAC OBSERVATIONS 

We make a continuous map of the central 0.77 square 
degrees of the Virgo cluster between M87 and M86 where 
many optical intracluste r light features have been iden- 
tified (iMihos et al.l l2005t iRudick et ahl l2009l see section 
i i5.2l and Figure [1]). We designed our observing program 
to reach the lowest possible surface brightness in the al- 
lotted time (100 hours; PID 60173). Our expected depth, 
based on the S'pzteer sensitivity estimator, is 0.0002 Mjy/sr 
at 3.6 fj-Tn. Discrepancies from this will be discussed below. 
An additional design constraint was that the map include 
an "off" cluster region for measuring the non-cluster back- 
ground level since we expect the diffuse intracluster light 
to be prevalent.^ 

Warm IRAC (jFazio et al.ll2004) simultaneously observes 
with two detectors at 3.6 & 4.5 /zm (channel 1 & 2 respec- 



tively). Each detector has a 5'.2 x 5. '2 field of view with 
1'.'2 pixels. Images are acquired in two adjacent fields, sep- 
arated by 6'.5. Our entire region is imaged in five layers of 
observations, each 10 min per pixel or 50 min per pixel to- 
tal over all layers. Inside of each layer we dither inside of a 
mapping pattern to cover the entire area while facilitating 
the removal of cosmic rays. We use a 100 sec frame time 
which gives a coverage of at least 30 overlapping frames 
per pixel. 3.6 & 4.5 /xm data were taken simultaneously 
and so have similar coverage, except for one array width 
on the north and south side of the final mosaic. Data 
were taken during two consecutive, two week long, IRAC 
campaigns PC15 & PC16 during February, 2010. Due to 
scheduling constraints, the AORs are spread out in time 
between these two campaigns, i.e. they are not taken con- 
tinuously in time. 

3. DATA REDUCTION 

Data reduction to reach our intended low surface bright- 
ness detection limits is a challenge because it is a factor 
of hundreds below the background level in the images and 
very little work has been done to characterize the behav- 
ior of the IRAC instrument at these levels. The following 
sections discuss our data processing steps with an eye to- 
wards the errors in each step. We start with a standard 
data reduction for both channels. We use the Spitzer Sci- 
ence Center (SSC) pipeline which generates corrected basic 
calibrated data files (CBCD; pipeline version S18.18) for 
our data reduction. These data have already been dark 
and flat-field corrected with standard calibration darks 
and flat-fields. A fiux conversion has already been applied 
based on A & K type standard stars placing the data in 
units of MJy/sr. 

3.1. Dark & Bias Correction 

Understanding the dark and bias contribution is impor- 
tant for LSB observations as it is a significant fraction of 
the total counts in each pixel. Calibration darks relevant 
to this dataset are taken as a set of 28 dithered 100s ex- 
posures within seven days of our observations. Standard 
darks can't be taken with IRAC because the shutter is 
not ever used. The darks are taken by looking at a re- 
gion of very low background near the north ecliptic pole. 
Those images are then median combined to remove sources 
which leaves an image of both the bias and dark current. 
Upon subtraction from the science data, both of these ef- 
fects are removed. Because the dark field does not have 
a completely zero celestial background, this dark removal 
technique will also remove some real celestial background. 
We account for this in our program design by intention- 
ally observing an "off" background location in the cluster 
which we will define to be of zero background, and tie the 
rest of the image to that (see H3.7I) . 

3.2. Flat Field Correction 

The accuracy of any LSB measurement is limited by 
fluctuations on the background level. A major contrib- 
utor to those fluctuations is the large-scale flat-fielding 
accuracy. The IRAC flat-field is derived from dithered 
observations, accumulated over 110 hours, of the zodiacal 
background as the best approximation to a source which 
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uniformly and completely fills the field of view. The pixel- 
wise accuracy of the flat-fields is very high, 0.2% & 0.1% 
at 3.6 & 4.5 ^m respectively. There is no variation with 
time in the flat field to the levels quoted. 

3.3. Illumination Correction 

There remains an illumination pattern in the science 
data in part due to latent images. Bright, saturated 
sources on prior images can fill traps in the arrays which 
decay away at a rate that is longer than the time between 
exposures thus leaving apparent flux at the location of 
that bright source in a previous image. Latent images can 
also result from the telescope slewing across a bright star 
which will leave a stripe of populated traps across the im- 
age. We removed these latents from the data by making 
a tailored illumination correction for each set of observa- 
tions. In each individual CBCD, we mask all the detected 
sources based on a SExtractor (jBertin fc ArnoutsI |1996[ ) 
segmentation map. SExtractor is run with a 1.5 sigma de- 
tection and analysis threshold and a minimum of 5 pixels 
for object detection. Then, the masked images are me- 
dian combined inside of each observing set to create the 
illumination correction. Median combining these masked 
images does an excellent job of creating an illumination 
correction because there are of order 100 dithered images 
per channel. Because we only want to remove the pattern 
in this illumination image, and not the background level, 
we first subtract the mean level from the illumination im- 
age, then we subtract the illumination pattern from the 
science data. 

3.4. Zodiacal Light Removal 

The dominant, non-instrumental background source in 
space-based infrared observations is the zodiacal light. 
This is light scattered from interplanetary dust and de- 
pends on both time of year and direction of observation. 
The zodiacal light is assumed to be uniform on the size 
scale of the IRAC field of view. A single value of the zo- 
diacal light per image is removed from the images using a 
model based on COBE DIRBE data (jKelsafi et al.lll998[ ) 
as listed in the headers. This model is expected to be ac- 
curate to 2% or better. Values for the zodiacal light are in 
the range of 0.104 to 0.114 MJy/sr at 3.6 ^m and 0.238 
to 0.323 MJy/sr at 4.5 /im. Zodiacal light values are listed 
in the headers in units of MJy/sr. However the conversion 
from data units to MJy/sr provided by the Spitzer Science 
Center is based on point sources. Because the zodiacal 
light is an extended source, this point source derived con- 
version is not appropriate and we apply a multiplicative, 
extended source aperture correction to the zodiacal light 
value (0.91 and 0.94 at 3.6 & 4.5 ^m respectively) before 
subtraction. 

3.5. First Frame Effect 

The IRAC arrays are known to have a relation between 
the bias level and the delay time between frames. This 
is referred to in the documentation as the "first frame ef- 
fect", even though it effects all frames. Here, delay time 
refers to the time between one exposure and the next, and 
not the exposure time itself. There are different delay 
times between our images which all have the same expo- 
sure time because of the dithering and tiling pattern we 



have chosen. Figure [5] shows this effect for our dataset. 
Specifically, 3.6 /xm data shows a logarithmic trend be- 
tween bias level and delay time. The mean level shown is 
measured from masked individual CBCD frames that are 
column pulldown corrected. Masks are derived from SEx- 
tractor segmentation images. We fit the first frame effect 
separately for PC15 and PC16 with a logarithmic function 
and apply this as a correction to the data. 

There is no measurable first frame effect in our 4.5 
/xm data (see Figured]). This is consistent with the known 
behavior of the arrays. 

3.6. Clock Time Effect 

After removal of the first frame effect, we notice a sec- 
ondary, undocumented effect between mean levels and 
clock time at 3.6 /zm. Figure |3] shows mean levels ver- 
sus frame number for each of the 25 AORs as well as the 
zodiacal light model pattern for the same frames. Spikes 
in this plot occur when a bright star or bright galaxy are 
in the frame. There are two correlations: (1) a ramping 
in mean level from the beginning of the AOR to the end, 
and (2) mean level variation from AOR to AOR. We see a 
large jump in background levels between campaign PC15 
data and campaign PC16 data, but also small level changes 
within the campaigns. 

The zodiacal light pattern is not likely to be responsible 
for the ramping effect. This is demonstrated in the lower 
portions of Figure [3] where the zodiacal light contribution 
is shown. It has a significantly different shape than the 
ramping effect. Neither the jump or the ramp are corre- 
lated with position on the sky or time since dark frames 
were taken. This rules out astrophysical background and 
build-up of flux on the chip as causes of the effect. 

It is likely that these effects are caused by the previous 
observation history of the chip. The ramping from the 
beginning to end of the AOR is likely a relaxation in the 
chip. As these AORs were not all taken sequentially in 
time, there is no way of knowing the precise observation 
history of the chips, and therefore no obvious correlations 
can be drawn with time since bright star observations or 
high background observations. The cause of the jump in 
background levels is unknown although potentially caused 
by zodiacal light model inaccuracies or previous observa- 
tion history as with the ramping. 

Because these effects can not be extragalactic in origin, 
we have chosen to remove them from the data. Each AOR 
trend is fit with an exponential function which is then sub- 
tracted from the data. We choose an exponential function 
as the appropriate model for a relaxation effect which fits 
the data well. 

At 4.5 /im there is the level change effect although not 
the ramping seen at 3.6 /xm (see Figure [3]). Again, we 
remove the level change from the data to put all AOR's 
on the same background level. 

3.7. The Mosaic & Background Determination 

After removal of the above discussed instrumen- 
tal and zodiacal light effects, we use MOPEX 
(|Makovoz fc Marleaull2005l) for the final astrometry align- 
ment and mosaicing, see Figure [TJ M86 is marked for 
reference. M87 is in the lower left corner, not imaged in 
our survey. 
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We determine a global background level from the north- 
ern edge of the mosaic as an ofF-cluster region so that we 
can have an absolute measure of ICL flux in cluster re- 
gions. This assumes we have taken care of all instrumental 
and local celestial effects, and the only remaining source of 
flux in the background regions is cluster light. This corner 
of both the optical and IRAC images appears free of low 
surface brightness cluster flux. Of course there still could 
be some cluster flux at this location, so any measurement 
of the diffuse cluster light from this data will be a lower 
limit. This spot has no saturated stars, is 30 arcminutes 
from M86, and 1.5 degrees from M87. Before measuring 
the background, we use SExtractor to make a segmenta- 
tion image. All sources are masked from the region. We 
also include only those pixels with a coverage of at least 
25 100s images. The background level determined in this 
corner is subtracted from the entire image to reach our 
final background level. 

3.8. The IRAC PSF 

Because our images are plagued by source confusion, we 
study the shape of the PSF to assure that there is not a 
lot of flux at large radii from all the objects which could 
overlap to mimic a LSB signal. We use a very bright, sat- 
urated star {K2mass = 5.5 mag.) in our mosaic to trace 
the outer wings of the PSF with high signal to noise, see 
Figure 131 This profile is made by first masking all objects 
in the frame, except the bright star, to twice their SEx- 
tractor measured isophotal radii. We then measure the 
median flux in radial bins out to ~ 300 arcseconds. For 
reference the dotted line shows an r~^-^ profile offset from 
the data for ease of viewing. Because the core of this star 
has been corrected for saturation, we can turn the radial 
profile into an encircled energy profile by integrating under 
the data curve. We find that 99.7% (99.6%) of the star's 
fiux is contained within 10 arcseconds of the central pixel 
at 3.6 fim (4.5 /xm). There is not a large amount of flux 
pushed to large radii by the PSF. 

4. DETECTION LIMITS 

In order to observe intracluster light features we bin our 
mosaiced data to achieve lower surface brightnesses. Bin- 
ning the data should reduce the noise in the image linearly 
with binning size. In Figure [5] we show a plot of noise ver- 
sus binning size for the data. We measured the surface 
brightness in nine (six for 4.5 /xm) regions of one pixel by 
one pixel across the SExtractor masked mosaic which are 
not near bright galaxies, stars, or optically detected ICL 
features. The noise is taken to be the standard deviation 
of those nine (six) regions. We then increase the size of 
the box incrementally, measuring the standard deviation 
at each step, until the box length is many hundreds of 
pixels. We are unable to measure this for a larger num- 
ber of regions because of these large final sizes and the 
available "empty" space in the image. Specifically the 4.5 
& 3.6/im images have different final shapes and so differ- 
ent numbers of blank regions. The actual length of the 
box is calculated from the square root of the number of 
unmasked pixels. For reference the solid line in Figure [S] 
shows a linear relation. The 'x' marks the sensitivity we 
expected to reach based on calculations with the Spitzer 
online estimation tool. The 3.6 & 4.5 iim mask shapes are 



calculated separately from each image to insure that the 
masking is consistently done to the same surface bright- 
ness threshold. The measured noise does not achieve the 
expected linear relation with binning length (see the Ap- 
pendix for further discussion of how to reach the lowest 
possible surface brightness limits). 

We believe that this technique for measuring the noise 
is conservative. Ideally we would like to measure the noise 
from many regions in close proximity to the plume in ques- 
tion which do not have any diffuse ICL in them. This is 
simply not possible as the diffuse intracluster light is ubiq- 
uitous and the area in the images limited. 

4.1. Small Scales 

There is a discrepancy at short binning length scales of 
just a few pixels. This discrepancy occurs because we are 
measuring noise on a mosaiced image so on small pixel 
scales few pixels), the noise is correlated and will not 
bin down appropriately. We have conflrmed that on a sin- 
gle CBCD image, the noise on small scales does behave 
linearly with binning length down to a resolution element. 

4.2. Medium Scales 

On scales from five to 30 arcseconds much of the ex- 
tra noise is due to sources in the image. This is demon- 
strated by changing the size of the masks. The first level 
of masking is just to use the SExtractor segmentation map 
as as mask (see ^3.31 for a description of parameters used). 
The resulting noise properties are shown with asterisks. 
This level of masking removes 16% of the pixels from the 
roughly 38.5 million pixels. Increasing the size of the 
masks to 1.5 (2.0) times the SExtractor determined ob- 
ject radii produces noise properties shown with a square 
(triangle) symbol. These increases in mask size increase 
the fraction of masked pixels to 62% and 75% respectively. 
Further increases in object mask radii only make inconse- 
quential changes to the noise properties, indicating that at 
twice the object radii is a good way to remove all of the 
flux from the galaxies. We also perform the test of using a 
mask that is the composite mask of 3.6 & 4.5 /im object 
shapes. This slightly changes the levels of the masks in this 
medium scale region in the same way that changing the in- 
dividual mask size affects noise measurement. Because we 
do see a change in noise level with changing mask size, we 
assume the discrepancy in this binning length regime is 
dominated by noise from the wings of galaxies which are 
improperly masked. 

Even after increasing the mask sizes there remains ex- 
tra sources of noise preventing us from achieving ultra-low 
surface brightness limits. There appears to be a floor to 
the noise at roughly 0.0005 MJy/sr/pixel at 3.6 /im and 
0.0008 MJy/sr/pixel at 4.5 /im. 

4.3. Large Scales 

We expect some of the large scale noise to be caused 
by the mapping pattern used for the observations. Our 
coverage map is not completely uniform, so on scales of 
roughly half a field of view there are differences in the to- 
tal exposure time and hence the total number of electrons 
detected. While this plot is shown in MJy/sr, we have 
performed the test of converting it to electrons using the 
coverage map. With this test we find an inconsequential 
change in the noise properties on large scales, certainly 
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nothing that brings it in line with the expected hnear re- 
lation. 

There are remaining sources of noise on these large 
scales, both instrumental and astronomical, which are very 
hard to disentangle. Uncertainties in the flat-fielding, and 
in our method of removing the first frame effect are two 
instrumental effects which are contributing to the noise 
on large scales. We expect the first frame effect to have 
a column-wise dependence that requires special calibra- 
tion data to measure. There is potentially also noise due 
to the ICL being blue in the IRAC bands while the zo- 
diacal light from which the fiats are made is red in the 
IRAC bands. Astrophysically, there is real structure in 
the zodiacal light and Galactic cirrus. There is also docu- 
mented diffuse intracluster light in the Virgo cluster itself 
and a small signal from the extragalactic background light 
which are both adding signal at low levels. Differentiating 
between these sources of error is no easy task (especially 
given the lack of a close-able shutter on IRAC). It is be- 
yond the scope of this paper to discuss the exact contribu- 
tions to the noise from these sources. Instead, we take this 
low surface brightness floor simply as an accurate measure 
of the noise on large scales. This floor level is very sim- 
ilar to the quoted extragalactic backgro und light level in 
publi shed IRAC work on that topic (See iKashlinskv et al.l 
l2005l also for a discussion of the errors in each of the 
potential instrumental and astronomical sources ). 

5. RESULTS & DISCUSSION 

Before discovering the noise floor, our program was de- 
signed to reach depths comparable to the LSB optical de- 
tections in this field with many resolution elements per 
plume. While we are limited by this fioor, there is still 
interesting science to be done with the IRAC data. A 
first-look through the images reveals ^ 100, 000 sources in 
both channels. 

5.1. A New Galaxy Cluster Behind Virgo 

Among the ^ 100, 000 detected sources, we have discov- 
ered a potential galaxy cluster behind the Virgo cluster at 
a redshift of z ~ 0.4. This candidate cluster jumps out as a 
grouping of galaxies in the IRAC data which is not notice- 
able in the optical imaging. We use SDSS catalog photo- 
metric redshifts to investigate this apparent over-density 
of sources in the IRAC images at 186.76731, -1-13.1463 
(J2000). Figure m shows the redshift distribution in a circle 
of radius 0.04 degrees centered on the candidate cluster. 
This distribution is compared to the mean redshift dis- 
tribution in nine similarly defined regions within the same 
area on the sky. The standard deviation on the comparison 
distributions is also shown. The apparent over-densities at 
z ~ 0.15 and 0.7 are less than two sigma. The peak at z 
~ 0.4 is highly significant at 11 sigma. This redshift peak 
is not detected in any of t h e published SPSS cluster cat- 
alogs (|Koester et al.l l2007t iMiller et al.l 120051: IWen et al.l 
120091 ). The only coincident objects listed in NED are the 
galaxies themselves and one ROSAT detected X-ray point 
source. This source requires further follow-up for confir- 
mation, including photo-z determination with the IRAC 
datapoints included, and spectroscopic confirmation. 

5.2. Intracluster Plumes 



Our goal in this paper is to detect and characterize the 
optically detected plumes in the region between M86 and 
M87. Our Spitzer images over lap with four optically de- 
tected regions identified bv iRudick et al.l ([2010) and shown 
in Figure lU These regions include three plumes(A, B, & 
D) and one region in the outer halo of M87 (El). Fluxes 
are measured by summing pixels inside of the four regions 
on our final reduced mosaic. As with the noise measure- 
ment, all objects on the final mosaic are masked to a radius 
of twice their detected radii (see i j4.2l ) . Surface brightness 
and noise measurements for these regions are shown in Ta- 
blCn At 3.6 fim, plume D is a five sigma, secure detection, 
while plumes A & B are not solid detections (1.9 and less 
than 1 sigma respectively). Region El is detected at nearly 
three sigma. At 4.5 /im, no regions are detected with high 
significance. 

Plume D is an inte resting source, likely Galactic in origin 
(jCortese et al.|[20Tol) . which will be discussed in a future 
paper. We focus here on plumes A & B and region El. 
Plumes A & B are not detected in the 100 fim IRAS im- 
ages. A non-detection implies that these plumes are extra- 
galactic, although it is still possible that they are Galactic 
features below the sensitivity of IRAS. 

5.2.1. Stellar Mass & Age 

We estimate both the stellar mass and age of the in- 
tracluster plumes by using our the individual facets of 
our optical and infrared data where they have the most 
power. There are three important pieces of information 
we can use. (1) B-V color puts a constraint on the age 
of the population. This is because the slope of the Wien 
side of the blackbody is more sensitive to age variations 
than the Rayleigh- Jeans. (2) IRAC flux upper limits put 
constraints on the upper limit of the stellar mass. The 
power of the infrared is in using the 3.6 /im flux to study 
mass. This is because most of the mass in an old popula- 
tion is emitting near the peak of the stellar distribution, 
whereas the young massive stars that emit heavily in the 
optical have faded away. This leads to large variations in 
the optical M/L ratio as a function of age, especially for 
old populations, while the infrared M/L is relatively con- 
stant over the same age range. Additionally, the infrared 
is a good place to study mass because it is not affected 
by dust extinction, as the optical is, and we do not expect 
hot dust emission to be affecting the IRAC fluxes due to 
the predicted old age of the population. Lastly, (3) Opti- 
cal fluxes put constraints on the lower limit to the mass. 
While the optical is not the best place to determine the 
stellar mass because of the reasons listed above, given the 
age of the stellar population from the colors, and assuming 
a dustless model, we can put a lower limit on the stellar 
mass. Any stellar population model that falls below our 
optical data points is ruled out by the optical detections. 

With these three pieces of information we constrain both 
age and stellar mass of the populations for plumes A & B 
starting with a color magnitude diagram of B - V vs our 
upper limits on 3.6 /im flux (FigurejTj). The 3.6 /im fluxes 
for plumes A & B are both considered upper limits to the 
actual fluxes because they are not more than three sigma 
above the noise. For the case of plume B, the flux mea- 
sured is less than one sigma, so we use the noise value 
itself as the observed flux density. For comparison, we 
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use a family of lBruzual &: CharlotI (|2003D models of a sub- 
solar metallicty ([Fe/H] = -0.3 3) simple stellar popul ations 
as elliptical galaxy templates. IWilliams et alj ()2007f ) work 
on HST ACS data of individual ICL stars in Virgo suggest 
they have sub-solar metallicity. The only variables in the 
models are age (1, 3, 5, 8, and 12 Gyr; our error bars do not 
warrant higher resolution) and mass(l x 10* to 4.5 x 10^ 
Mq). We choose to use this elliptical template since ICL 
sources have been shown to be optically red, both from 
observations (|Krick et al. 2006; Rudick e t al.l 20Tol) and 
from theoretical work (jWillman et al.ii2004l : Murante et al.l 
120041) . 

We first examine age from the B-V color. We expand 
the work of iRudick et al.l (j2010i) by showing the full range 
of model ages consistent with their measured B-V col- 
ors for the plumes A & B. Plume A is consistent with a 
population older than 3Gyr, while plume B is consistent 
with a population older than 5Gyr. Age determinations 
from the optical colors measured are also affected by the 
metallicity of the plumes. In Figure |7| we show the direc- 
tion that the mass curves would move if the metallicity 
were increased by a factor of 2.5 to the solar value. If 
the metallicity is larger than our assumed value, then the 
plumes are younger than we predict, and viceversa. 

We next examine an upper limit to the stellar mass from 
the 3.6 fim fluxes. We quote here a range of masses 
which fit inside the error bars of the ICL plumes, they 
represent a conservative upper limit to the mass of the 
plumes. The upper limit to the mass of plume A is 
1.1 X 10^ - 4.5 X 10^ M©. The upper hmit to the mass 
of plume B is 3. x 10*^ - 1.5 x 10^ Mg. We have also 
checked that the iMarastonI (|2005l ) models give consistant 
estimates. 

Lower limits to the mass are inferred from to com- 
bination of the age model with the optical fluxes. We 
again look at the simple stellar population model described 
above. We choose the lowest age that fits the B-V color 
(3Gyr for Plume A, and 5Gyr for plume B), and deter- 
mine which masses are ruled out by the detections in the 
optical bands which gives a lower limit to the masses of 
the plumes. Using a model without dust here will give the 
lowest lower limit because assuming dust will only increase 
the inherent optical fluxes and thus increase the mass esti- 
mates. The lower limit to the mass of plume A is 5.5 x 10* 
and B is 2.1 x 10*. In conclusion, assuming the highest 
upper limit and the lowest lower limit we find the masses 
of plume A and B to be between 5.5 x 10* — 4.5 x 10^ and 
2.1 X 10* — 1.5 X 10^ M0respectively. These are the first 
mass measurements of intracluster light plumes. 

The masses we find for plumes A & B are consis- 
tent w ith those found in the simulations of iRudick et al.l 
(|2009[ ) . Other theoretical studies of ICL production do not 
discus s plume/stream prop erties (He nriqucs & Thomas. 
20101: iMurante et all 120071 : ISommer-Larsen et all 120051 
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For comparison we calculate the total stellar mass 
in Virgo from the sp ectroscopic luminosity function of 
iRines &: Gelled (|2008f ). We integrate under the Schechter 
function fit to their luminosity function from Mr = — 1 1 to 
Mr = —23. This gives us a total luminosity of 2.9 x 10^^ 
Lq. If we assume a mass to light ratio for this cluster of 
500 (iSchindler et all [19991) . we find a total cluster mass 



of 1.5 X 10^^ Mg. Th is is in the same range as found by 
ISchindler et"all (|l999t) based on X -ray and optical obser- 
vations who find a total mass of both the M87 and M49 
groups of 3 X 10" ± 0.6 Mq. 

In total, if plumes A & B are extragalactic, and if they 
were the only intracluster light in the Virgo cluster they 
would represent at most a couple percent of the total stel- 
lar mass in the cluster. Our Spitzer mosaic covers roughly 
half of the area around M87 and half of the area around 
M86/M84, so if the density of plumes/streams is constant 
in the region around the center of the cluster we could 
imagine that at most a few percent of the cluster mass is 
in ICL streams at this particular point in the Virgo clus- 
ter's evolution. Therefore, at any given time in the clus- 
ter's history, the plumes do not account for a large fraction 
of the cluster mass; but we expect plume generation and 
destruction to be ongoing in the cluster. 

Simulations suggest that these plumes las t for less than 
one Gyr, or roughly 1-2 crossing times (jRudick et al.l 
|2009| ). Assuming the cluster has been forming for roughly 
ten Gyr, and assuming a constant rate of plume forma- 
tion, we naively expect to find a few tens of percent of 
the total stellar mass be formed in plumes over the life- 
time of this cluster. This is consistent with total ICL 
fractions as measured from the diffuse ICL component of 
10-40% of the total cluster light fKrick fc Bernsteinl l2007l: 
iGonzalez et~al]|2005l : iFeldmeier et al. 20041 ). This suggests 
that bright plumes can account, order of magnitude, for 
all of the diffuse ICL as measured. We then don't need 
to rely on ram pressure stripping (RPS), or other gas dy- 
namical effects for input to the diffuse ICL. Although these 
processes do occur, they might not be the dominant input 
mechanism to the ICL. This could explain why intraclus- 
ter light appears red, without any evidence for in situ star 
formation which you might expect to see if RPS were the 
dominant input mechanism. While we can't rule out the 
other enrichment mechanisms, our data is consistent with 
gravitational processes which create bright plumes as the 
main enrichment mechanism to the ICL. 

5.3. M87 Halo 

We combine our infrared data on the halo of M87 to un- 
precedented radii with both optical and infrared archival 
data of the center of M87. The infrared data will allow us 
to understand if there is a dust component in M87 which 
may be effecting the optical estimates of age. From this 
age gradient information we hope to learn about the po- 
tential formation mechanisms of this giant elliptical. 

To learn about the stellar populations in the outer halo 
of M87, we examine the optical and infrared color pro- 
file to large radii in Figure jS] The outer halo of M87 is 
defined by region El, shown in Figure [H 160 kpc (2008 
arcseconds) from the center of M87. This is four times the 
radius analyzed by a recent multi-band optical study by 
iLiu et all (|2005l ) (even though the profiles in their plots do 
extend to comparable radii) . Opti cal measurements o f the 
same outer ha l o regi on comes from lMihos et al.1 (|2005l ) and 
IRudick et al] (l2010D . For comparison, both infrared and 
optical ph otometry of the central region come s from the 
literature (|Shi et al.ll2007l: IZeilinger et al.lll993D . Because 
M87 has a central jet, we specifically choose photometry 
from the literature that carefully removes the jet compo- 
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nent. 

For reference, the color-color plot includes a Grasil 
model evolution ary track for an elliptical galaxy 
(iSilva et alJ 119981) for ages ranging from formation to 13 
Gyr. These models are different from standard stellar pop- 
ulation models because they include dust in three environ- 
ments; interstellar HI clouds heated by the general inter- 
stellar radiation field, star forming molecular clouds and 
HII regions, and dusty envelopes of dying stars. These 
models are intended to be elliptical templates. They have 
an initial infall of gas with an exponential decay timescale 
of O.lGyr followed by passive evolution. We choose to use 
models with dust in them because dust is important in pre- 
dicting the infrared colors of a younger population. Note 
that these models show infrared colors becoming bluer 
with age, which is opposite the color trend of the optical 
data. When we discuss trends below we will use "blue" in 
the sense of the optical data. 

We use our data to study the existence of dust in clus- 
ters. Dust grains should be present in the ICM because 
we know that metals, which are often depleted onto grains, 
escape galaxies by ram pressure stripping, supernova ex- 
plosions, and tidal interactions It is not clear how long dust 
is able to survive in the ICM due to the destructive effects 
of sp uttering off of the hot ICM ()Yamada fc Kitavamal 
l2005t ). The direct detection of dust in the far-infrared has 
so far not been possible due to the large background fluc- 
tuations and limited telescope resolution, making it dif- 
ficult to d etermine if si gnal is from th e ICM or galaxies 
(iBai et all [2 007: Kitay ama et all l2009t iMontier fc GiardI 
120051 : Giard et al.l 1200^ However, there is substantial 
evidence for intracluster dust. Many low-redshift cluster 
galaxies exhibit significant amounts of extraplanar 8 m 
emission (e.g. NGC 4501 & NGC 4522) which is thought to 
be the result of polycyclic aromatic hydrocarbons (PAHs) 
begin pulled out of the galaxy by ram pressur e strip - 
ping. Further evidence comes from iGirardi et all (|1992[ ). 
who find that asymmetric redshift distributions in nearby 
groups implie s extin c tion c o nsistent with dust o bscuration. 
Additionally, Mao2 (IT995h : iMcGee fc Ba"ioghl (^2010) find 
that high-redshift objects behind clusters show large vi- 
sual extinctions. 

To get a broad idea of the presence of dust, we compare 
age predictions from the optical and infrared data. Age 
predictions from our infrared colors are consistent with 
age predictions from the optical data within 2a which im- 
plies broadly that dust is not effecting the optical colors. 
Dust can cause extinction at optical wavelengths. With a 
low surface brightness measurement, there is no good way 
of ruling out dust extinction without infrared data which 
can easily penetrate the dust. 3.6 & 4.5 ^m observations 
provide this without being effected by hot dust emission. 

The second thing to notice is that any interpretation of 
a color change in the profile from the inner region to the 
outer halo at 160 kpc is a two sigma measurement. Both a 
flat color profile, and a slightly blueing optical profile are 
consistent with our measurements within two sigma. In 
other words, our data do not place strong constraints on 
the color profile of the halo of M87. 

Optical work on the color profile at large radii has 
been interpreted as having a blue trend attributed to 
younger ages with radius and/or a potential metal- 



licity gradienllZeilinger et al.l ()1993[ ): iLiu et al.l (|2005l ): 
iRudick et a\\i201(lt) . We note that the detection of a blue 
trend in the optical work at large radii is also a low sigma 
measurement. 

A flat or blue optical gradient in the outer halo is 
also consistent with simulations. Specifically, ICL simu- 
lations form the ICL stars at redshifts greater t han 1.5 
fMur ante et al.ll200 4': Som mer-Larsen et al.ll2005D . which 
would therefore make them red at low redshift, caus- 
ing the outer halos of galaxies like M87 to also be red 
lik e their centers. Nume rical simulations of two clusters 
in ISommer-Larsen et al.l ((2005) predict a fiat to slightly 
blue-ward profile in optical colors from a radius of lOkpc 
to IMpc from the center of the brightest cluster galaxy 
(BCG). They find a decrease in metallicity with radius 
causing the blue- ward trend. According to this work then, 
differentiation between a blue or flat optical gradient in the 
outer halo would give us information on the metallicity of 
the ICL population. 

Since this color profile extends to the largest radii yet 
discussed, we are likely sampling the ICL population. The 
relatively old stellar ages in both the inner and outer pro- 
file to the limits of our observations suggests that the halo 
of M87 has not been infused with a young population in 
the last few Gyr. A young ICL population could come 
from either in situ formation triggered by interactions with 
other galaxies or the cluster potential. Or it could come 
from the stripping of already young stars away from an 
infalling galaxy. The red colors of the outer halo rule out 
these cases for this massive galaxy in the deep potential 
well of center of the Virgo cluster. 

We hesitate to say that the flat color gradient in the 
outer halo is a byproduct of monolithic collapse. Because 
the surface brightness profile shows a distinct, outer ICL 
component |Liu ct al. 2005), it is likely that the outer com- 
ponent has a different formation mechanism than the inner 
deVaucouleurs profile. However, our findings are consis- 
tent with the outer profile being composed of stars simi- 
larly old as those in the center of the galaxy. 

6. CONCLUSION 

In summary, we describe Spitzer IRAC imaging of the 
central 0.77 square degrees of the Virgo cluster with spe- 
cific attention paid to the low surface brightness features 
in this region. This is the first intracluster light study in 
the infrared. We reach a limiting SB of 3.6 x 10^ — 4) 
MJy/sr for the largest regions. We don't detect the in- 
tracluster light plumes as seen in the optical. From up- 
per limits on their brightnesses we present an optical and 
infrared color magnitude diagram and the first ever mea- 
surement of the mass of an intracluster light plume. We 
overlay that plot with a family of model elliptical galaxies 
which allows us to put an upper limit on the mass and age 
of the plumes. Then optical colors are used to constrain 
the lower limit to the masses. We do this for two plumes 
finding masses in the range of5.5x 10^ — 4.5x 10^ and 
2.1 X 10* — 1.5 X 10^ M0for the two plumes for which we 
have coverage. These are consistent with streams formed 
in simulations. In total, these types of plumes, assumed 
constant formation rate and summing over the lifetime of 
the cluster, are massive enough to account for the entire 
diffuse ICL population. 
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Secondly we look at the color profile of MS 7. Age pre- 
dictions from our infrared colors are consistent with age 
predictions from the optical data which implies that dust 
is not effecting the optical colors. Comparing our infrared 
data on the outer halo of M87 with literature values for the 
inner region and optical profile, we find that the M87 color 
profile at large radius is either flat or becomes slightly more 
optically blue with radius. A blue profile would indicate 
younger or lower metallicity stars at larger radii. Although 
a flat color profile is strictly consistent with monolithic col- 
lapse models, M87 is known to have a diffuse outer ICL 
component with a different radial profile, implying a dif- 
ferent formation mechanism compared to the rest of M87. 

Lastly, we comment on the intricacies of the IRAC in- 
strument, and the best methods to observe LSB features. 
LSB work with IRAC requires careful attention to details 
of the instrument and of the astrophysical background. 
Due to instrumental effects, large scale binning of the data 
does not allow the observer to reach deeper surface bright- 



ness detection limits as expected. Increased exposure time 
does afford the deepest possible surface brightnesses, as 
described in the appendix. 
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ysis Center/California Institute of Technology, funded by 
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California Institute of Technology under a contract with 
NASA. This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and 
Space Administration. 
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Fig. 1. — Top: The full 3.6 mosaic, including intraclustcr and M87 halo regions as used in the text. Bottom: A zoom in on plume D 
where the low surface brightness signal is visible in the Dl region. 
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Fig. 2. — First Frame Correction: All ~ 3000 images at both 3.6 (squares+triangles) and 4.5 fim (asterisks + pluses) are shown in this 
plot of delay time vs. mean background level. Mean levels have been arbitrarily scaled to show both channels on the same plot. The data 
from the two campaigns have different mean background levels as shown here with the different symbols within each channel. For the 3.6 
/im data the mean value in each bin of delay time is shown with the white asterisks which is then fitted with a logarithmic function shown 
by the solid lines. 
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Fig. 3. — Relaxation effect for 3.6 fim (left plot) and 4.5 ^m (right plot): Each of the 25 AOR's is shown in this plot of frame number 
within the AOR vs. mean background level (arbitrary units) as solidly squiggly lines. AORs within campaign pcl5 are shown in black, and 
those in campaign pcl6 are shown in blue. Zodiacal light model levels are shown as smooth lines in the bottom third of the plot, arbitrarily 
scaled to fit in the frame. Spikes in mean level occur when a bright star or galaxy is on the frame. Our cycling dither pattern exposes some 
objects on multiple consecutive frames leading to spikes wider than one frame number. 
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Fig. 4. Top; IRAC PSF to 300 arcscconds. 3.6 & 4.5 /iin are shown in solid and dashed hnes respectively. The dotted hne shows an 
^-2.5 cm-yQ for reference. Bottom: Encircled energy of the above IRAC PSF, here showing only the last 1% of the stars flux. The IRAC PSF 
is well behaved and does not push a lot of flux to large radii. 
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Fig. 5. — Noise properties of the mosaic; 3.6 fj,m on top, 4.5 //m on bottom. For reference, the right hand y-axis is in units of nw/m^/sr 
used in extragalactic background light studies. Asterisk symbols axe for simple object masking, squares and triangles show masking to 1.5 
and 2.0 times the object radii. The 'x' marks the sensitivity we expected to rea<;h based on the Spitzer online performance estimation tool. 
The solid line shows a linear relation for reference. While object masking does bring the noise levels down, we still do not rea<;h the levels 
expected from binning implying further noise sources creating a noise floor at ~ 0.0005&0.0008 MJy/sr/pixel for 3.6 & 4.5 /nm respectively. 
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Fig. 6. — Redshift distribution of a galaxy cluster behind the Virgo cluster. The solid line shows the redshift distribution of all SDSS 
detected galaxies with photometric redshifts (and errors on those rcdshifts below 0.2) within 0.01 degrees of the center of the new cluster. The 
dashed line shows the mean redshift distribution in the same size regions near to this new cluster. The error bars are one standard deviation 
on the mean of those background distributions. 
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Fig. 7. — Color magnitude diagram of plumes A & B for stellar mass and age determination. Points are the data for both plumes including 
one sigma error bars or upper limits. Solid lines show constant mass tracks as a function of age. Masses are listed at the oldest point on 
each line. Ages are shown for the 1.0E9 mass track. An arrow indicates the direction that the mass tracks would move if the metallicity were 
increased to the solar value. 



14 



Krick et al. 




Fig. 8. — Color-color plot of the M87 halo. The points indicate the colors of the inner region and outer halo of M87. The line is a Grasil 
model of elliptical galaxy evolution as a function of age. The outer halo is consistent with having the same color as the inner region of the 
galaxy. 
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Table 1 
Plume Properties 



Plume ID 


Signal 


Unmasked Length 


Real Length 


la 


noise 


# of sigma 




MJy/sr 


AB mag. 


arcsec 


arcsec 


MJy/sr 


AB mag. 




A 


0.00067 


14.28 


322 


668 


0.00036 


0.65 


1.9 


B 


-0.00039 


15.78 * 


203 


410 


0.00045 


0.76 


< 1 


Dl 


0.00256 


15.16 


111 


228 


0.00048 


0.21 


5.3 


D2 


0.00155 


15.05 


156 


308 


0.00046 


0.34 


3.4 


El 


0.00129 


14.68 


186 


401 


0.00046 


0.41 


2.8 


A 


0.00065 


14.32 


347 


668 


0.00052 


1.2 


1.3 


B 


0.00043 


15.56 * 


218 


410 


0.00055 


0.78 


< 1 


Dl 


0.00128 


15.91 


118 


228 


0.00068 


0.65 


1.9 


D2 


0.00030 


16.84 


166 


308 


0.00062 


1.2 


< 1 


El 


0.00111 


14.84 


197 


401 


0.00059 


0.65 


1.9 



* These values for the signal are actually the noise values used since this is a less than one sigma detection 

Note. — Top half for 3.6/im, bottom for 4.5jUm. The unmasked length for the two channels is not equivalent 
because we measure the unmasked area in each channel individually. Real length refers to the size scale of the 
plume as defined in Rudick et al. 2010 
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APPENDIX 

HOW TO DETECT LOW SURFACE BRIGHTNESS FEATURES WITH IRAC? 

We have shown above that there is a "floor" in reaching low surface brightnesses when binning together many pixels 
because the noise does not drop as expected with binning length. We now investigate if lower surface brightnesses can be 
achieved by increasing exposure time., i.e. does t he noise decrease w ith the square root of exposure time as expected? To 
look at this we have used IRAC dark field data (jKrick et al.ll2009t ) as a unique dataset that goes very deep in the same 
region of the sky mainly devoid of background. We don't want any signal to be in the background (low zodiacal light, 
low Galactic diffuse emission, no intracluster light) so we can remove as many astrophysical sources of noise as possible. 

Using all the warm mission dark calibration data to date, we have made a mosaic of 300 frames, each with 100s exposure 
time. Mopex was used to make this mosaic from the column pulldown corrected CBCDs. We mask all sources with a 
SExtractor segmentation image. We look at all the unmasked pixels with the same coverage levels, ranging from 12 - 
300 images. The noise on the distribution of background values within each coverage level is the standard deviation of 
the Gaussian fit to that distribution. Each distribution has > 750 pixels in it. For comparison we have done this same 
analysis on the dark field mosaics from the first year of the cryogenic mission. 

These plots show that IRAC noise does decrease roughly as expected with exposure time. The slight deviation at larger 
exposure times is likely caused by the first frame effect and by source wings. The dark field data has not been given the 
special treatment of deriving it's own first frame effect as discussed for the Virgo dataset above. Also, the sources have 
been only simply masked with the SExtractor segmentation map, and so wings of sources will still add noise to the images 
which affects the measurement at lower surface brightness. 




Fig. A9. — Noise as a function of increasing exposure time; 3.6 fim on left, 4.5 fim on the riglit. Stars and solid lines are the warm data 
and fitted square root function. Squares and dashed line are the same for the cryogenic mission. 



